Non-paraxial Airy beams 
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We report the propagation dynamics of Airy light beams under non-paraxial conditions. It is studied using 
the general approach which deals with Fourier expansion of the beam. We show the transformation of the 
beam from the Airy form as the paraxiality parameter decreases. The role of evanescent waves in non-paraxial 
regime is discussed. © 2009 Optical Society of America 
OCIS codes: 260.0260, 350.5500. 



During some last years it was reported [1] about 
exciting experimental realization of the Airy light 
beams. There were indeed confirmed the important non- 
diffracting and accelerating properties of these beams. 
From the theoretical point of view the Airy beams are 
described by the diffraction equation [2] 



du 1 d 2 u 
~d~z + 2k Uz 2 



= 0, 



(1) 



which is valid only for paraxial waves. 

In the current Letter we study the properties of non- 
paraxial Airy beams, which are the exact solutions of 
the Maxwell equations. We will perform the theoretical 
analysis of propagation of non-paraxial beams and this 
can offer the new problems for experimentalists in the 
field. 

Let us consider the light beam propagating along z- 
direction in isotropic medium with dielectric permittiv- 
ity e and magnetic permeability \i. We suppose that in 
the initial plane z = the field is distributed according 
to the Airy function as Ai(x/xo) exp(ax/xo), where a is 
a decay factor limiting the beam energy at x < 0, x$ 
is an arbitrary transverse scale. Then this field changes 
during propagation according to the Maxwell equations. 

If a = 0, we get to the ordinary Airy light beam, which 
is non-diffractive as it was discovered in the pioneering 
work [3]. However, such a beam should possess the infi- 
nite energy to be created and, therefore, cannot be real- 
ized experimentally. Parameter a confines the beam to be 
generated, so that it is not more non-diffractive, but can 
be experimentally generated. Note that there are other 
possibilities to confine the beam [4] . 

It is well-known that an arbitrary beam distribution 
in the plane z = can be constructed using the plane 
waves. Each plane wave propagates as it follows from the 
Maxwell equations: 

H(x, z, k x ) = (a(k x )e y - c 2 (k x )^^j e ^+^ z , 



/ e x k\ 

B(x,z,k x )= f c 2 {k x )e y + ci(fc a ) v j 



, (2) 



where the amplitudes ci and c 2 (in general, complex 
numbers) describe polarization of partial plane waves. 
The dependencies on the transverse wavenumber k x in 
the amplitudes c\ and c 2 describe the profile of the 
beam. If ci = (c 2 = 0) then wave is TE (TM) po- 
larized. Transverse k x and longitudinal k z wavenum- 
bers are connected with the dispersion equation, so that 
k z (k x ) = ^/k^e^i — k 2 , where ko = ui/c is the wavenum- 
ber in vacuum. Also, the wavevector is of the form 
k = k x s x -\- k z Q z . 

In further the dimcnsionlcss parameters are intro- 
duced: transverse wavenumber q = k x xo, wavenumber 
in vacuum x = ^o^Oj transverse i = x/xq and longitu- 
dinal z = z/koXy spatial ranges. 

To obtain the Airy light beam in initial plane z = 
we assume that the amplitudes are equal to [2] 



ci,2(<?) = 77--4i,2e" 

Z7T 



^i(q 3 —3a 2 q—ia 3 ) /3 



(3) 



where Ai and A 2 are complex constants to provide an 
arbitrary beam polarization. 

Now the Fourier transform (the superposition of the 
plane waves incident at different angles) is for the propa- 
gation of a non-paraxial Airy beam. For example, for the 
electric field of the beam one reads 

1 f°° ( Ai 
E (£) S ) = 7T A ? e y + — (-? e * + C(g)ex) 

27T J -co \ £ X 



xe 



iqx+ixC(q)z -aq 2 i(q 3 -3a 2 q-ia 3 ) /3 



dq, 



(4) 



where ((q) = \Jx 2 e[i - q 2 . 

This beam is the superposition of the TE and TM Airy 
beams. Both of them are described in the same manner. 
Therefore, we can investigate only TE Airy beam. TM 
beam is characterized by another amplitude coefficient 
A. TE Airy beam has only y electric field component 
and follows from Eq. when Ai = 0. 

Paraxial Airy beams, which are usual for experiments 
[2,5-7], can be obtained, when the spectrum of propagat- 
ing waves is wide. It corresponds to the high-frequency 
radiation, which is described by the paraxiality param- 
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Fig. 1. The distribution of normalized intensity of the 
Airy beam under paraxial condition (x = 100) calculated 
using Eq. (J4j) . The decay parameter a = 0.1. Refractive 
index is n = 1. 



eter x » !• Then one can use the approximate ex- 



pression a/x 2 £M- <7 2 ~ Xv 7 ^ 2 X V£M' 
solution ((4]) transforms to the paraxial one 



and the exact 



Ey(x, Z) 



2tt 



xe 



-aq 2 e i(q 3 -3a 2 q-ia 3 )/3^ 



(5) 



where n = y^e/I is the refractive index. The integral ([5]) 
can be expressed using the closed-form expressions [2] as 

E v (x, z) = A 2 e mx2? -Ai [x - (5'/2) 2 + iaz'] 
xcxp (ax - az' 2 /2 - i(z' s /12 - a 2 z' /2 - xz'/2)) , (6) 

where z' = z/n. The equivalence of Airy beam de- 
scriptions given by Eq. ^ for large parameter x an d 
paraxial formulae ([5]) and ([6]) can be directly verified. 
Fig. [1] shows the distribution of normalized beam inten- 
sity \E y (x, z)\ 2 /\A2\ 2 calculated from Eq. Q in paraxial 
case (for great number x) ■ One can see the typical struc- 
ture of Airy beam and its main property - acceleration - 
and compare them with previously reported results ob- 
tained on the basis of Eq. ([6]) [2,5]. 

As the paraxiality parameter x decreases, the beam 
structure becomes more and more different from that in 
Fig. [TJ When coordinate z is small, the non-paraxiality 
weakly influences (see Fig.[2ja)). The acceleration of the 
non-paraxial Airy beam is the same as that of the parax- 
ial one, because it is determined by the amplitude spec- 
trum ([3|). Only after z = 2 the beam starts to diverge 
rapidly. At first, the side lobes disappear, while the main 
lobe is still strong. After z ~ 4 the non-paraxial Airy 
beam dissipates. The discussed Airy beam has the non- 
paraxiality parameter x comparable with the unity. The 
properties of the paraxial beams quickly restore, if x 
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Fig. 2. The distribution of normalized intensity of the 
non-paraxial Airy beam calculated from Eq. (j4|) at (a) 
X = 5, (b) x — 2, (c) x = 1) (d) X = 0.5. Parameters: 
a = 0.1 and n = 1. 
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Fig. 3. The initial (at z = 0) transverse distribution of 
normalized intensity of the non-paraxial Airy beam at 
(a) x = 1 an d (b) x = 0.1 carried by propagating waves, 
evanescent waves, and both propagating and evanescent 
waves (total beam). Parameters: a = 0.1 and n = 1. 



increases. For example, for x — 10 the spatial range, 
where the acceleration is excellently seen, significantly 
raises (approximately by two times) . When x further de- 
creases, the beam's acceleration is getting less and less 
pronounced (sec Fig.[2](b), (c), and (d)). Now the beams 
propagate almost strongly ahead, while the side lobes 
rapidly attenuate. 

The explanation of the behavior of the non-paraxial 
Airy beams can be provided with dividing the spectrum 
of the plane waves forming the beam into two parts. The 
first one describes propagating waves and corresponds to 
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Fig. 4. The distribution of normalized intensity of the 
non-paraxial Airy beam at % = 0.1 corresponding to (a) 
propagating waves, (b) evanescent waves, (c) total beam. 
Parameters: a — 0.1 and n = 1. 



the waves with real longitudinal wavenumber k z , there- 
fore, —x^feji < q < Xy/zJi- The other qs (q < —Xy/eji 
and q > —X\f^V) are f° r the evanescent waves, which de- 
cay during propagation because of imaginary wavenum- 
ber k z - When parameter x is great, the spectrum of the 
propagating waves is wide, while evanescent waves can 
be neglected. That is why one uses the paraxial approx- 
imation for the field of the Airy beam ©. When x is 
comparable with the unity, both propagating and evanes- 
cent waves should be taken into account. In Fig. [3] we 
show the field distributions of propagating and evanes- 
cent beams, when the total beam energy corresponds to 
the Airy beam. For x — 10 (not shown) the profile of the 
propagating beam coincides with that of the total (both 
propagating and evanescent) intensity distribution. The 
maximal value in the profile of the evanescent waves is 
less than 10 -10 , therefore, the evanescent waves can be 
certainly neglected. For x = 1 the contributions of the 
evanescent and propagating waves are comparable. The 
energy of the propagating waves is concentrated in the 
main lobe, that is why it propagates at long distance (see 
Fig. [2] (c)). The energy of the evanescent waves is local- 
ized in the side lobes. This is the reason, why the side 
lobes rapidly decay. In Fig.[3](b) the beam with very nar- 
row spectrum of the propagating waves is shown. We can 
observe very low energy of the propagating waves. This 
beam can be approximately considered as the evanes- 
cent Airy beam. Note that the sum of the energies of 
the propagating and evanescent waves is not the total 
energy of the beam, because the last includes also the 
interference term (it arises due to the interference of the 
propagating and evanescent waves). 

In Fig.|4]the propagation of the almost evanescent Airy 



beam is shown. The low intense propagating beam prop- 
agates without spreading at the shown distance. How- 
ever, it yields very low input to the total intensity. In 
spite of the suppressed intensity of the evanescent beam 
its main lobe propagates at long distance without sig- 
nificant change of the intensity. It should be noted that 
the realistic distance of beam's propagation decreases 
with diminution of paraxiality parameter x-i because 
z = X x oZ- 

As for the energy flux of the TE and TM Airy beams, 
it has x and z non-zero components. More interesting sit- 
uation arises for a superposition of these partial beams, 
which realizes a specially polarized Airy beam. Then the 
Poynting vector possesses also y component. However, 
this case require separate investigation. Similar study 
for the Bcssel beams was earlier performed in Ref. [8]. 

In conclusion, we have proposed the non-paraxial Airy 
beam as the exact solution of the Maxwell equations. 
We have theoretically analyzed the propagation of this 
beam and interpreted its properties using those of prop- 
agating and evanescent waves. We have proved the non- 
paraxiality parameter x to be large enough (at least 
X > 10) to provide the paraxial beam and its proper- 
ties (acceleration). If x is small, the evanescent waves 
prevail and the beam rapidly decays. The beam with 
small x can be regarded as the evanescent Airy beam. 
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